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Abstract 

We review some of the theory relevant to passage times of one-dimensional Levy 
processes out of bounded regions, highlighting results that are useful in physical 
phenomena modelled by heavy-tailed Levy flights. The process is hypothesised to 
describe the motion of a particle on the line, starting at 0, and exiting either a hxed 
interval [—r, r], r > 0, or a time-dependent, expanding, set of intervals of the form 
rf^], r > 0, K > 0. Asymptotic behaviour of the exit time may be as r 1 0 or 
as r ^ oo, but particular emphasis is placed herein on “small time” approximations, 
corresponding to exits from or transmissions through thin strips. Applications occur 
for example in the transmission of photons through moderately doped thin or thick 
wafers by means of “photon recycling”, and in atmospheric radiation modelling. 

Keywords: Levy flights, Levy processes, alpha-Stable processes, photon transmission, 
thin-wafer semiconductors, passage times. 

1 Introduction 

Levy flights are increasingly being proposed as models for physical phenomena in a va¬ 
riety of applications areas. Convincing empirical evidence for their existence is reported 
for example in experiments or observations on the transmission of photons through scat¬ 
tering media such as clouds, especially relevant to climatological studies (e.g., i) and 
the transport of “holes” in semiconductors (e.g., m)- The latter two applications, in 
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particular, are concerned with the one-dimensional motion of a particle or object in a 
bounded domain - a “slab” in the case of [3], a “wafer” in the case of by means 
of jumps whose magnitude has a heavy tailed distribution, specihcally, a power law tail. 
The corresponding random walk - the cumulative sum of the jump sizes - has an infinite 
variance but is in the domain of attraction of a stable law. Since stable (Levy) processes 
in continuous time can be obtained as limits of interpolated Levy flights, we can expect 
the particles’ motions to be well described by a stable process in one (or more) dimen¬ 
sions. In many contexts these provide better descriptions of the process than a Brownian 
motion model. 

In the present paper we review some of the theory relevant to passage times of one¬ 
dimensional Levy processes out of bounded regions with a view to highlighting results 
useful in these kind of modelling contexts. The results will be formulated in terms of 
general random walks and Levy processes taking values in the real line M, but special case 
reductions to (one-dimensional) stable processes and Levy flights will be made along the 
way. 

A process Xt, t > 0, starting at 0, is hypothesised to describe the motion of a particle 
on the line, starting at 0, and exiting either a hxed interval [—r, r] or a time-dependent, 
expanding, set of intervals of the form [—Here k > 0 is a parameter to be 
specihed and r > 0 is a variable which in the physical context (and in the case k = 0) 
would specify the thickness of the slab or wafer. Introducing the time dependence allows 
for the slab thickness to increase in dimension at an algebraic rate in time, starting from 
zero thickness. Of course in the physical context this would represent a limiting case 
where the particle’s escape is restricted in some way. 

Exits of Xt from the regions just described are represented by passage times of the 
form 

Tnir) = inf{f > 0 : \Xt\ > rt^}, r > 0. (1.1) 

Expressed in space-time coordinates, Ti^{r) is the time at which the point {t,Xt) first falls 
outside the rectangular (when k = 0) or curvilinear (when k > 0) region dehned by 

Rk '■= t > 0, a: G M, |a:| < rt^}, (1.2) 

having started at (0, 0). When k = 0 this corresponds to exit from a hnite slab of width 
2r. 

There are two possible asymptotics of interest. Most common in applications is to 
take r oo, a “large region”, hence “large time”, approximation. In [3], for example, the 
application is to photon free path distribution in atmospheric radiation modelling, where 
the heavy tailed distribution of path lengths arises with the passage of light through cloud 
types where rapid aerosol density fluctuations result in free paths varying from less than 1 
meter to many kilometres. “(Levy-stable) distributions are therefore well suited to model 
photon free paths globally” m)- 

At the other end of the scale, in m, short wavelength radiation is shone onto doped 
crystal wafers of thicknesses of the order of 50-350 pm, creating an electron-hole pair. 
On a nanosecond time scale, the hole moves randomly with the thermal velocity (this 
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motion well described by Brownian motion) until recombining with an electron. Emitted 
photons are absorbed, thereby generating a new hole, from which a new photon is emitted. 
The resultant “photon recycling” propagates through the crystal and its location, when 
projected on the r/-axis, dehnes a 1-dimensional random walk. But in this instance the 
jump distribution is better described by a power law tail decay having inhnite variance. 
After a large number of jumps a stable distribution is to be expected for the particle’s 
location, and this is indeed observed experimentally. “The power law decay of the hole 
concentration ... is steep enough at short distances ... to £t the data for the thin sample, 
and at the same time slow enough at large distances ... to account for the data for thick 
samples” (M, p.20). Here, for thin samples, concern is with the short time behaviour of 
the approximating stable process and the asymptotic of interest for T^^ir) is as r j, 0. 

In the next sections we consider both “thin strip” (short time) and large time be¬ 
haviour. First we need some preparatory material in the next section. 


2 Functional Laws for Random Walks 


“Short Time” only makes sense for continuous time processes. How to formulate these 
in a way that’s relevant to the physical situation? We can start with a random walk S'„ 
comprised of i.i.d jumps i = 1, 2,...; thus, 

n 

Sn = J2^^, n = l,2,..., (2.1) 

i=l 


with Sfj = 0. Each has the distribution of a generic rv which is assumed to have a 
heavy tail satisfying 

F(|^| > x) ~ as X —)■ oo, (2.2) 


as well as the tail “balance” condition 

P(^ > x) 

lim —--- exists and equals c, 0 < c < 1. 


(2.3) 


In fl2.2p . CK is a parameter in (0,2). Conditions fl2.2p and 02.31) imply that the random 
walk is in the large-time domain of attraction of a stable random variable, S^°‘\ of index 
a, which means that there are centering and norming sequences An and Bn, such that 


Sn-Ar, 

Bn 


D 


^(a). 


as n — )■ OO. 


(2.4) 


Here and below denotes convergence in distribution. When 0 < a < 1 in 02.2p . the 

An can be taken as 0, and when 1 < a < 2, F|^| is hnite and we can set An = nE{^), 
which we can further take as 0 by recentering. Further, when 02.21) and 02.3p hold, the 
Bn satisfy Bn ~ as n —?• oo for some constant a > 0 which can be taken as 1 by 

rescalingpj 

^When a = 1 tlie centering is logarithmic. More general norming sequences, where is multiplied 
by a slowly varying function, rather than a constant, can be catered for in a more general setup. But the 
situations we consider here suffice to illustrate the main points. 
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as n —>■ oo. 


(2.5) 


Assume then that we have the convergence 

D ^ n(a) 

^ ’ 

This extends to a “functional” version as follows. For each f in a hnite interval [0,T] 
let Sn(t) be the interpolated process dehned by “joining the dots” in the graph {j,Sj), 
j = 0,1,2,... ,n. Then we have, in an appropriate sense, the convergence 

t Si \ as n ^ CX), (2.6) 

where (5'^“^), f > 0, is a stable process of index a. 

We can dehne discrete-time passage times by 

T^(r) = min{n =1,2,...: \Sn\ > rn^}, r > 0, (2.7) 

and in view of the functional law in fl2.6p . they will converge as n —)■ cxo, after rescaling, 

to the corresponding passage times for dehned by 

T^^\r) = inf{f > 0 : > r^}, r > 0. (2.8) 

Having adopted a Levy hight (S'„)„=i^ 2 ,... or a stable process (*S'}“^)t>o as a physical model, 
properties of the time taken for the process, starting at 0 at time 0, to traverse and exit 
the regions dehned in fll.2p . can then be derived. In the continuous time process we can 
analyse small time behaviour, as f j, 0, and this translates for the passage time to small 
“slab width” behaviour, i.e., as r j, 0. Large time behaviour, as r —)■ oo, can similarly be 
considered. For the discrete time model, only large n behaviour can be considered. 

There are two ways to think about the modelling procedure. One, is that we model the 
physical process by the discrete time random walk Sn and think of the stable process as 
a continuous time approximation to it. Alternatively, we could simply model the physical 
process directly with the continuous time stable process. Much of the time the difference 
will just be one of interpretation because in practice the continuous time model has to be 
htted via a discretised version in any case. 

3 Functional Laws for Levy Processes 

Rather than the discrete time Levy flight, we may adopt as a physical model a general 
Levy process {Xt)t>o and concern ourselves directly with distributions and behaviour as 
r j, 0 or as r —>■ cxo of the passage times Tf^{r) dehned in fll.ip . Just as for the discrete 
time random walk, a real-valued Levy process (W)t>o ni^y be in the domain of attraction 
of a stable law either as f j, 0 or as t —)■ oo. Conditions on the Levy tail characterising 
this are available in [12] (the canonical triplets of the rescaled underlying processes are 
required to converge to that of the stable process in an obvious way). If these conditions 
obtain, there is a functional law analogous to fl2.6p which gives that a rescaled version of 
Xt converges to a stable process (S'^' ^)t>o either as f j, 0 or as t —)• cx). 
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4 Basic Stable Process Properties 

Stable processes, t > 0, are pure jump Levy processes whose distributions at each 
hxed time t > 0 are stable with index a in (0,2). The distributions are absolutely 
continuous; thus, the distributions of 5^“^ at each hxed time t > 0 possess densities. 
Moments of 5'^“^ of orders less than a are hnite but moments of orders greater than 
or equal to a are inhnite. Particularly attractive for physical modelling is that stable 
processes have the scaling property: for each A > 0, has the same distribution, 

as a process, as In particular, for each A > 0, has the same distribution 

as so in a sense which can be made precise, 5^“^ is of order A^/" both for small 

and for large A. Note also that the scale invariance implies a form of self-similarity of the 
process. 

A Levy process is completely determined by its canonical triplet ( 7 , 0 -^, IT), where 
7 G M is a shift constant, cr^ > 0 allows for a Brownian component (for the stable laws, 
= 0), and n is a measure on M which determines in a certain sense the rate and size 
at which the jumps of the process occur (for a Brownian motion, there are no jumps, and 
n = 0 .) n is in turn determined by its positive and negative tails, which for the stable 
laws are of the form 

^ > 0 ’ ( 4 - 1 ) 

Jb Jb 

where c± > 0 are constants with c = c+ -f c_ > 0 , and a G ( 0 , 2 ) is the index of the stable 
law. We refer to 

n(a;) = n^(a;) + n"(a:) = ^, x > 0 , (4.2) 

simply as the tail of the measnre. An important point is that 11 possesses a singularity at 
X = 0 of order a; this corresponds to a concentration of inhnitely many small jumps w.p.lH 
in any hnite interval. One can see the possible relevance of this property for modelling in 
the thin wafer sitnation. 

For further properties of stable processes we refer to [13] and [T5] . 


5 Passage Times out of Rectangular or Curvilinear 
Regions 

In this section we consider the passage times Tf^{r) of a general Levy process dehned in 
fll.ip . but including special reference to the stable processes. To begin with we note a 
distinct difference between the cases k = 0 and k > 0. When k = 0 and r > 0 the region 
is a rectangular (slab or wafer) of width 2r not including (0,0). Bnt when k > 0 the 
curvilinear region shrinks to 0 at small times. In the second case, there is the possibility 
that the continuous time process might jump out of the region instantaneously, even for 
a hxed value of r > 0 . At another extreme, there is the possibility that the process never 

^“w.p.l”= “with probability 1”, or almost surely 


5 





gets out of the region at all; in this case we formally set Tf^{r) = oo. The latter is never 
the case when k = 0 but can occur when k > 0 since the width of the intervals then 
expands indehnitely. We need to keep these cases in mind when analysing the passage 
times. In general they hinge on the relative sizes of the two parameters n and a, as we 
now discuss. 

Instantaneous Exit 

To rule out the possibility of an instantaneous exit we have the following equivalences, 
true for any Levy process X, which can be obtained by elementary calculations. For any 
r > 0 , 


P{Ti^{r) > 0) > 0 P{Ti^{r) > 0) = 1 liminf < r, w.p.l. (5.1) 

40 P 

Instantaneous exit is thus ruled out when the righthand relation in fl5.1l) holds. We can 
get information on this from [2], Theorem 2.1. From it we deduce that for k < 1/2, 
lim^o = 0 w.p.l for any Levy process, so the righthand relation in fl5.ip always 

holds in this case; while for n > 1 / 2 , lim^o \Xt\/P = 0 w.p.l if and only if (iff) n(x'^)da: 
is hnite. Referring to fl4.2p . we can conclude: Instantaneous exit is ruled out for the stable 
process 5*^") if 0 < n < 1/a. The interpretation is that as long as the singularity in fl at 
0 is not too steep {a < 1 /k), takes a nontrivial time to escape from R^. 


Exit in Finite Time 

When TK(r) = 0, exit obviously occurs in hnite time, but when > 0 there is 

the possibility that = oo, which is to say that w.p.l X never escapes from the 

region in hnite time. This situation can be analysed as follows. It is convenient to dehne 
X/ := supo<;s<i \Xs\/s'^, t > 0, with X^ := sup^^g \^t\/P- Then for f > 0 

iUr) >t}C {X < r} C {T.(r) > t}. (5.2) 


In the Appendix Section 0 some equivalences are proved for P{TX) < oo) = 1) and it is 
shown that this holds ih P(|X/ > rP for some f > 0) = 1. Sufficient for this is that 


lim sup 

40 


P 


> r, or lim sup 

t—^OO 


P 


> r, w.p.l. 


(5.3) 


fl5.3|) describes escape of X^ from inhnitely often w.p.l either at small times or at large 
times. But it is easier to deal with than the equivalences in flQ.ip of the Appendix. 

Theorem 2.1 of can be used again for information on the lefthand condition in (15.311 . 
For K < 1/2, it never holds, while for k > 1/2, we can deduce that limsup^^o \^t\/P > 0 
w.p.l ih limsup^io \^t\/P = oo w.p.l ih n(a:'^)da: is inhnite. For the stable process 
the latter is the case ih 0 < k < l/a. 

For the righthand condition in (15.3p . we refer to |7], which gives limsup^^o^ \Xt\/P = 
oo w.p.l except when k > 1 or 1/2 < k < l,P|Xi| < oo and EXi = 0. In these latter 
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cases, limsupj_^j^ \Xt\lt'^ = oo w.p.l iff n(a:'')dx is infinite. Specialising to the stable 
process the integral diverges iff k < 1/a. Putting together both sides of (15.dh . we 
conclude: 5'^“^ escapes from in finite time w.p.l for any k > 0 and 0 < a < 2. 

Moments of the Passage Times 

Restricting ourselves at hrst to the discrete time situation, we can deduce from Prop. 
2.2 of Doney and Griffin (2003) that ET^{r) is hnite for any r > 0 when 0 < k < 1/2. 
The same property carries over to the continuous time version. 


6 Asymptotics 

The results in the previous section apply for each hxed r > 0. In a physical situation we 
expect the process will only approximately behave as a stable process. The approximation 
may be close (and desired to be close) at small times (corresponding to small slab widths r 
for the passage times) or it may be close at large times (corresponding to large widths r). 
If a Levy flight (discrete time heavy tailed random walk) is assumed for the motion, fl2.4p 
and fl2.6p describe approximation by a continuous time stable process. If the underlying 
model is taken to be a general (continuous time) Levy process, it may also be well described 
by a stable process for small or large time. 

In either situation, because of the functional laws, the passage times in continuous 
time (dehned by fll.llL or in discrete time (dehned by fl2.7l) ). converge after appropriate 
rescaling to the corresponding passage times of the approximating stable process either 
as r 0 or as r —)■ cxo, corresponding to small or large time approximations. These have 
the properties listed in the previous section, and others related to the distributions of the 
passage times. 

6.1 Relative Stability of Passage Times 

“Relative Stability” refers to the convergence in distribution of Ti^{r), either as r 0 
or as r —)■ cxo, to a nonzero degenerate random variable. Thus, convergence is in fact 
in probability to a hnite nonzero constant. This is a special case of the convergence in 

p 

distribution of Tf^{r). Specihcally, we have relative stability when Tfir)/C{r) —)■ 1, for 
some nonstochastic function G(r) > 0, where the convergence is in probability as r 0 or 
as r —)■ oo. 

This kind of behaviour occurs iff the underlying process X is “relatively stable” in the 
sense that Xt converges in probability to a hnite nonzero constant, either as t 0 or as 
t —)■ oo. This property has been thoroughly analysed and conditions for it to occur are 
listed in P and HD]. They hold in particular when X has bounded variation with positive 
drift (for r 0), or when X has hnite, positive mean (for r —)■ oo). In these cases the 
norming function C{r) described in the previous paragraph can be taken as a constant 
multiple of 0 < k < 1. When k = 0, corresponding to exit from a slab of width 

2r, this means that Tfir) behaves approximately linearly as r 0 or as r —)■ oo. 
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6.2 Passage Times Distributions 

Returning to the random walk model defined in fl2.1l) . a relevant large-time functional law 
for the passage time is in Doney and Mailer (2005). It tells us that when fl2.6p holds 


T^jr) D 1 

j-i/ Yi/u ’ 


as r ^ oo, 


( 6 . 1 ) 


where u = 1/a — k and n < 1/a. In flO.ip . Y has the distribution of supo<s<t 

The restriction k < 1/a ensures that Y is finite w.p.l. Similar relations hold for the 

continuous time ^^(r) defined in fll.ip . 


7 Overview 

The following overview may help to distinguish the various processes discussed in the 
paper. We may have as model: 

• a discrete time random walk in the domain of attraction of a non-normal 

stable law as u —)■ oo (i.e., the random walk is a “Levy flight”) 

• which converges as n —)■ cxo, after rescaling, to a non-normal stable process ^^)t>o 
(via the functional law in fl 2 . 6 p h so that 

• the discrete time passage time in fl2.7p converges to the stable process passage time 
in fl 2 . 8 p (again, via the functional law). 

Alternatively, we may have as model: 

• a continuous time Levy process (W)t>o in the domain of attraction of a stable law 
either as t 0 or as t —>■ oo 

• which converges either as f 0 or as t —)■ oo to a stable process (5'^“^), t > 0 , (via a 
corresponding functional law), so that 

• the continuous time passage time in fll.ip converges to the stable process passage 
time in fl2.8p (via the functional law). 

Alternatively, again, the continuous time model may be the stable process ((S^^^), itself. 

Having constructed passage times T^{r), T^ir) or Ti“^(r) from one of these processes, 
we can study the limiting behaviour as r —)■ oo (for any of T^(r), Ti.^{r) or Ti“^(r)) or as 
r I 0 (for Tc(r) and Ti"^(r)). 


8 Summary 

There is a vast literature on domains of attraction both for continuous and discrete time 
processes, and for stable distributions and processes themselves. We have discussed only 
the 2 -sided exit problem as this seems to be relatively newly considered in the applications 
literature. One-sided exits, i.e., from a semi-infinite domain, have been more commonly 














studied and applied, for example in the fluctuation theory of random walks and Levy 
processes. Still, many gaps remain to be hlled. 


9 Appendix: Notes and Further References 

For exit in hnite time we have the following equivalences. Recall the dehnition := 
suPo<s<i \Xs\/s'^, t > 0, with := sup^^o Then 

P{T^{r) < cx)) = 1 lim P{T^{r) > t) = 0 lim P{X^ < r) = 0 

t^OO t^OO 

< r) = 0 P(sup \Xt\/P < r) = 0 -^(1X^1 > rP for some f > 0) = 1. 

t>o 

(9.1) 


Some information on the rate of divergence of ET^{r) as r —)■ cx) is in |1] and [5]. 
Acknowledgements This research was partially supported by ARC Grant DP1092502. 
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